Multifractal zero mode for disordered graphene 
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Off-diagonal disorder with random hopping between the sublattices of a bipartite lattice is de- 
scribed by a Hamiltonian which has chiral (sub-lattice) symmetry. The energy spectrum is symmetric 
around E = and for odd total number of lattice sites an isolated zero mode always exists, which 
coincides with the mobility edge of an Anderson transition in two dimensions(2D). In the chiral 
orthogonal symmetry class BDI we compute the fractal dimension D2 of the zero mode for graphene 
samples with edges. In the absence of disorder D2 = 1, which corresponds to a one-dimensional edge 
states, while for strong disorder D2 decays towards and the zero mode becomes localized. The 
similarities and differences between zero modes in the honeycomb and the square bipartite lattices 
are pointed out. 

PACS numbers: 72.80.Vp,71.30.-l-h,72.20.Ee,73.22.-f 



I. INTRODUCTION 



The interest in zero energy modes extends to more 
than a decade and to all areas of physics [1|. It was 
recently revived owing to the emergence of graphene [2| 
and topological insulators [3j. The zero modes are be- 
lieved to signal non-trivial lattice topology and are con- 
nected to excitations of fractional charge [1|. They also 
exist in the spectrum of a bipartite lattice (two inter- 
connected sublattices) if the total number of lattice sites 
is odd [J]. In the presence of off-diagonal disorder which 
connects the two sublattices, in the so-called chiral dis- 
ordered systems, a zero mode is a critical point (mobility 
edge) of the Anderson metal-insulator transition. The 
Dirac point of graphene is a stable zero mode in two 
dimensions {2D) which becomes a critical point for ap- 
propriate off-diagonal disorder, and zero modes with half- 
integer charge appear in superconducting Bogolyubov-de 
Gennes systems as Majorana fermions 5]. 

Graphene has attracted immense interest in the last 
decade. In this novel material the electrons can be de- 
scribed by Dirac equation and its band structure displays 
two valleys (act as pseudo spin states) related by time- 
reversal [3|. The Fermi surface has two non-equivalent 
Dirac points so graphene is a semimetal or a zero gap 
semiconductor, at long length scales it has a linear en- 
ergy dispersion with a density of states which goes to zero 
at E = Oy. On the other hand, disorder causes Ander- 
son localization for all states in 2D disordered systems, 
at least in the limit of infinite size[6|. This is expected 
also for graphene since the honeycomb lattice has a co- 
ordination number even smaller than the square lattice. 

Our interest here is the so-called off-diagonal disorder, 
which appears in the random hopping among the sites 
of two interconnected sublattices A, B (mirror images of 
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one another), of a bipartite 2D lattice such as the hon- 
eycomb lattice of graphene. For graphene off-diagonal 
disorder between nearest neighbours (n.n.) can appear 
by removing lattice bonds at random, it can simulate 
ripples, etc. It is expected (as diagonal disorder) to mix 
graphene's two-valleys and if correlated to break time- 
reversal symmetry by decoupling valleys, etc. The off- 
diagonal disorder (random hopping between A and B) is 
different from usual diagonal disorder (random site ener- 
gies) , in the sense that it respects the discrete chiral (sub- 
lattice) symmetry[7|. The chiral disordered systems were 
classified by three random matrix ensembles (the chiral 
universality classes BDI, AIII, CII)t8|. A zero mode can 
exist a.t E = 0[9| where the chiral system undergoes a 
metal-insulator transition in 2D driven by topologically 
induced Anderson localization[10]. This type of localiza- 
tion is shown via non-perturbative field-theoretic tech- 
niques and is controlled by vortex-like excitations of the 
sigma model, via a mechanism similar to the Berezinskii- 
Kosterlitz-Thouless transition[ll| . 

The chiral symmetry of the Hamiltonian H, with ran- 
dom hopping among the A, B sublattices, manifests itself 
in the energy spectrum which is composed of random 
energy pairs {Ej,—Ej), where j labels the eigenvalues 
of the same sign. For an odd total number of lattice 
sites the spectrum is symmetric around E = and al- 
lows zero modes which cannot match, as required by chi- 
ral symmetry, any of the positive or negative Ej^s. In 
fact, the E = mode matches to itself. A zero mode 
is a critical state which emerges among localized states 
for topological reasons, the localization length (and the 
density of states) diverges at ii^ = 12|, |l3| . In 2D off- 
diagonal disorder can localize all non-zero states, the only 
exception is the zero modefl^. We shall show that off- 
diagonal disorder can localize even the zero mode when 
it becomes strong [10[. In the original treatment of chiral 
systems [l2| no traces of weak localization (quantum in- 
terference) was found to all orders in perturbation theory. 
We show that the critical zero mode displays multifrac- 
tal fluctuations [15| described by a fractal dimension D2, 



which tends to zero as the disorder becomes very large. 
In other words, for strong ofF-diagonal disorder the zero 
mode state becomes locahzed in agreement with [1(| . 

Our purpose is to study the niuhifractal properties of 
the zero mode in the absence of magnetic field or spin- 
orbit coupling, that is for a Hamiltonian H with real 
random hopping among A and B which belongs to the 
chiral orthogonal class BDI. This H has a mobility edge 
in 2D which is a zero mode. In the absence of chiral 
symmetry in 2D a mobility edge can exist only in two 
circumstances, if the usual orthogonal universality class 
is replaced by the symplectic with spin-orbit coupling, or 
by the unitary class with added magnetic field. In the 
presence of disorder two mobility edges usually appear as 
critical points of an Anderson transition. For a chiral 2D 
system they merge into one at the zero mode whose origin 
is topological [9j. Our study has an advantage over 31?, 
where mobility edges usually exist, since they are more 
difficult to handle. For the chiral class BDI in 2D we 
compute the fractal dimension -D2 16[ of the zero mode 
state. This is done in order to see the fate of the Dirac 
point of graphene in the presence of off-diagonal disor- 
der. One might also ask the question about the difference 
of the obtained D2 between a honeycomb and a square 
lattice, since both lattices with appropriate off-diagonal 
disorder belong to the same chiral class BDI. Another 
question addressed in this paper is whether the critical 
zero mode localizes for strong disorder and the nature of 
its localization. 



II. THE 2D HAMILTONIAN WITH 
OFF-DIAGONAL DISORDER 



The honeycomb lattice consists of two interconnected 
triangular sublattices A and B. If the lattice has 
zigzag(zz) edges zero modes can appear located at the zz 
boundary of one sublattice[17]. A spinless electron with 
off-diagonal disorder is described by the single particle 
nearest neighbour tight-binding Hamiltonian (belongs to 
the chiral orthogonal class BDI): 



H 



r6Ai=l,2,3 



ialbr^ 



(1) 



Or, ^r+si annihilates (creates) an electron in A or i? sub- 
lattice. The A, B sublattices are interconnected by ran- 
dom hopping ir,i, r G A sublattice and Si,i = 1,2,3 con- 
nects to its three nearest neighbours (n.n.) which belong 
to the B sublattice. The nearest neighbour hopping t^ ^ 
are chosen real uncorrected random variables distributed 
by the law: 
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which guarantees positive t only (its mean < t >^ 0). 
This is equivalent to a box distribution in the range 
[—W/2,+W/2] for the logarithm ln(i). Our choice of 
n.n. off-diagonal disorder (its strength is denoted by W) 
respects the chiral (sublattice) symmetry, while for pure 
graphene tr^i = t for every r, i {W — O). A local sym- 
metric gauge model of off-diagonal disorder (also chiral) 
has zero mean < t >= 0, t not ln(i) is chosen, e.g. from 
a box distribution of width W , which gives a single value 
of off-diagonal disorder (no other energy scale exists in 
the Hamiltonian and for the mean zero hopping the lat- 
tice vanishes). In the gauge model all ensemble aver- 
ages are invariant under sign changes of the wave func- 
tion and the single particle averages are site-diagonal jl2| . 
The Hamiltonian of Eq.(l) for both models of n.n. off- 
diagonal disorder (logarithmic and zero mean) has no 
spin-dependence, it is invariant under time-reversal (t is 
real), and it respects chiral symmetry (belongs to BDI). 

A bipartite lattice consists of the two sublattices A and 
B, in the A-B basis the Hamiltonian is: 



H = 



H 



HAi 




AB 



(3) 



Hab contains the random hopping between A and B. 
The chiral symmetry can be expressed in the anti- 
H, (J3] = with the pseudospin 

For an energy E and corre- 



commutation relation 

[10 
matrix (J^ — 

sponding state ^ = 






H^ = i?*, via the above 



commutation relation an energy ~E and a state cts^ = 

*4 \ 
-T also exists, Hias'i') = -as.H'i! = -Ela^,^). 

-^B J 

The presence of chiral symmetry implies pairs of ener- 
gies E, —E with wave functions \1/, (73 \E', respectively. The 
spectrum is symmetric spectrum around E = and for 
an odd number of lattice sites a state with zero energy 



^A 






*A^0 

*i3 = 



(4) 



and its amplitude has support only on the A sublattice 
being zero on the B sublattice. 

The Hamiltonian H has at least Na — Nb independent 
zero modes [71 , Na is the number of sites for sublattice A 
and Nb the number of sites for sublattice B, with Na > 
Nb- The total number of lattice sites N = Na + Nb 
is odd and since the number of energies is also odd the 
chiral symmetry leads to an i? = mode. We choose 
Na = Nb + 1 for one zero mode to exist, which can be 
easily obtainedy, |l6| by solving: 



H^ ^0^ 



H 



H 



AB 



AB 








= 0^ 



Bab's B = 
HIb'^a - 
(5) 



The '^A are obtained from the homogeneous system 
H^g'i'A = of Nb hnear equations defined on the B 
sublattice, for the amphtudes on the A sublattice. The 
system has Na — Ng + 1 unknowns (the additional equa- 
tion is provided by the normahzation) [15| . 

The critical zero mode state of Eq.(l) is multifractal, 
it has infinite self-similar detail which roughly repeats it- 
self. In the renormalization group sense multifractality 
implies the presence of infinitely many relevant operators. 
The critical wave functions at the disorder induced An- 
derson transition are multifractals 15] . For weak disorder 
the periodic Bloch states become chaotic and for strong 
disorder (above a critical disorder for dimensions higher 
than two) become localized in a region of the order of 
the localization length ^. A finite ^ for a localized state 
can be estimated by the participation ratio (outside the 
localization region denoted by ^ the states usually de- 
cay exponentially). At the zero mode the localization 
length diverges logarithmically and precisely at -E = 
takes any value which depends strongly on the choice 
of boundaries [4;]. At the Anderson transition a critical 
wave function lies at the mobility edge which separates 
extended (chaotic) from localized states. The multifrac- 
tal states in the border between chaotic and localized 
states have filamentary structure and are described by 
multifractal dimensions [15]. In the chiral problem stud- 
ied, however, no chaotic (diffusive) states exist and the 
critical zero mode emerges among localized 2D states 
only. 

The multifractals are inhomogeneous fractals described 
by many (a whole spectrum) of fractal dimensions Dq, q e 
[—00,00], obtained from scaling the distributions for a 
physical quantity rather than its mean (ordinary frac- 
tals are characterized by a single fractal dimension D via 
scaling of the mean). Alternatively, the multifractals de- 
scribe a singularity spectrum f{a), it is related to D' s 
by an appropriate Legendre transform. In our case the 
irregularly distributed probability amplitude l'0i,m|^ at 
the lattice site l,m of the normahzed (^; ^ |V'i,m|^ = 1) 
zero mode state is collected independently of /, m. In 
the presence of chiral symmetry a state with energy E: 
ipi^m, I, ITT- — 1,2,... along the sites of the bipartite lat- 
tice, at the mirror image energy —E corresponds to the 
state: (— l)'^™^/';,™. The JV'/,™]^ is not self-averaging and 
can be described by the spectrum of dimensions obtained 
from scaling of all the moments of the wave function. In 
order to describe l'0;,m|^ of the zero mode we introduce 
multifractal measures (the moments of the wave func- 
tion), in other words one must compute the spectrum of 
dimensions 



for g = 2. The D2 (known as correlation dimension) 
describes scaling of the inverse participation ratio 



D„ 



1~ q L 



uni — , q e 



InL 



-00, 00 



(6) 



IPR^Y.l'f'^- 



(7) 



with respect to the linear length L of the square sample. 
In the rest we focus only on the fractal dimension D2 



for large L. The I PR gives information for the degree 
of localization (it is the inverse for the number of lattice 
sites participating in the state). The D2 is conveniently 
obtained by scaling the averaged < I PR > , or the typical 
exp < \n{IPR) >, vs. the linear size L, D2 is the slope 
in the corresponding log- log plot. For a metal a chaotic 
(extended) state has D2 — d, the space dimension, for a 
localized state D2 -^ (insulator), while for a multifrac- 
tal critical state (at the metal-insulator transition) D2 
lies in between these limits. The multifractal dimensions 
Dq for the wave function are closely related to scaling 
for the moments of the local density of states, e.g. for 
the zero mode they can represent correlation functions 
for the local density of states. 



III. THE ZERO MODE FOR THE HONEYCOMB 
LATTICE 



We have solved the linear equations of Eq.(5) for the 
honeycomb lattice which consists of Nc horizontal zig- 
zag chains (the width of the lattice is ^V^ in units 
of a = \/3ac-C: o.c-c is the distance between carbon 
atoms) and Ngc the number of sites for each chain (the 
lattice length is L = ^°^~"'" in units of a) . The lattice has 
N = NcNsc sites, N is odd if Nc and Nsc are odd. We 
have also taken NcV^ ^ Ngc — 1 which corresponds to 
almost square L x L graphene samples without bound- 
ary conditions (be). The chosen samples have zz edges 
in the two parallel boundaries and in the other two have 
armchair edges. The amplitudes of the zero mode wave 
function are obtained by building a statistical ensem- 
ble of random Hamiltonians for every value of the off- 
diagonal disorder W. For W — the lattice symmetry 
with zz edges guarantees the presence of zero mode edge 
states [17]. These zero mode states appear due to the lat- 
tice topology (the Hilbert space structure), for diagonal 
disorder see [18[. Thus, the zero mode is an edge state 
having amplitude at the zz boundary and fractal dimen- 
sion D2 = 1. For non-zero off-diagonal disorder W which 
respects the chiral symmetry the fractal dimension D2 
diminishes. For strong W the zero mode state localizes 
along the zz edges. 

In Fig.l we plot the distribution of the I PR for the 
zero mode state. In Fig. 2 the mean < IPR > and the 
geometric mean exp(< \n{IPR) >) are plotted in a log- 
log plot vs. the linear size L. From the slopes one can 
obtain the fractal dimension D2 , for zero disorder D2 = 1 
and for strong disorder (large W) it decreases to zero 
(Z?2 — >■ 0). For graphene the localization of the zero 




IV. COMPARISON WITH THE SQUARE 
LATTICE 



The above analysis applies to any bipartite lattice (also 
for the square) with appropriate off-diagonal disorder 
which respects chiral symmetry. This is true as long 
as appropriate boundary conditions are taken and the 
size allows presence of the zero mode[9'], e.g. vanishes by 
taking periodic boundary conditions. We repeated our 
computations for square lattice N = Lx L samples, with 
N odd for the E = mode to appear. In Fig. 3 the prob- 
ability distribution for the IPR is shown. In Fig. 4 the 
log-log plots of the mean and the typical from the slopes 
allow to compute the fractal dimension D2 for the square 
lattice. 



FIG. 1: The flow of the probability distribution P{ln{IPR)) 
obtained from half a million graphene squares of linear lengths 
about 1/ = 8, 10, 11, 13, 15, 16 with logarithmic off-diagonal 
disorder of strength W — 1. 



mode wave function occurs on the edge. For finite values 
of W^ ^ the fractal dimension D2 takes values between 
zero and one. 
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FIG. 2: The log-log plot for the scaling of the mean < IPR > 
and the typical exp < \n{IPR) > with graphene samples of 
linear size L. The data are taken for half a million realiza- 
tions with logarithmic off-diagonal disorder W = 1, and the 
straight lines are the fits which give the fractal dimensions 
D2. 



FIG. 3: The fiow of the probability distribution for the 
\n{IPR), taken from half a million realisations of square lat- 
tices of lengths L = 25, 41, 49, 61, 69, 81, 101, with off-diagonal 
disorder W = 1. 



In Fig. 5 are shown the results of D2 for both bipartite 
lattices, graphene and square. In the absence of disorder 
{W — 0) for graphene the zz edges contribute to a zero 
edge mode with dimension D2 = 1 while for the square 
D2 = 2. For low disorder, e.g. W = 0.1, D2 is close 
to one for graphene and two for the square, respectively. 
For higher disorder W the D2 take values between one 
and zero for graphene and two to zero for the square. 
For strong off-diagonal disorder the wave function be- 
comes localised along a one-dimensional path, e.g on the 
zz edges for graphene samples, and localized in the bulk 
for the square[16|. For strong disorder {W — 5) the frac- 
tal dimension D2 approaches zero in both cases. 

Our conclusions concern D2 for graphene and the 
square lattice in the presence of off-diagonal disorder 
both with chiral symmetry. The D2 obtained from scal- 
ing the average IPR is expected to be more sensitive to 
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FIG. 4: The log of the mean < IPR > and the log of the 
typical exp(< ln(7P_R) >), vs. the log of the linear size L. 
The off-diagonal disorder is W — 1 for half a million realiza- 
tions. The straight lines are the fits which give the fractal 
dimensions D2. 



rare events than scaling the typical value which turns out 
to be more reliable in this case. The two 132 's, derived 
form scaling the mean and the typical, are comparable 
with each other with rather small uncertainty (Fig. 5). A 
rather large difference is seen for disorder higher than 
W ~ 1. One could denote this disorder as a point of a 
freezing transition (difference between averaged and typ- 
ical values) [19|. The main difference between graphene 
and square occurs for zero disorder, where for graphene 
the zero mode state has a fractal dimension D2 = 1, 
which corresponds to a wave function localised along the 
zz edge (edge state), while for the square D2 — 2 for an 
extended state. For finite values of disorder the 132 's in 
graphene take values between one (the wave function is 
localised on the edge of the sample) and zero. For strong 
disorder both systems have D2 — 5> 0. 



V. DISCUSSION-CONCLUSIONS 



In this paper we have studied the multifractal prop- 
erties of the zero modes which appear in chiral orthog- 
onal systems (class BDI) for an odd number of lattice 
sites. The zero mode is simply a zero energy state of 
the real symmetric Hamiltonian with off-diagonal disor- 
der which respects chiral symmetry and occurs exactly at 
the critical energy (mobility edge) of an Anderson metal- 
insulator transition in 2D chiral disordered systems. This 
transition has only insulating and critical (E — 0) phase, 
it has no metallic phase. The fractal dimension D2 for 
the critical state at i<J = is computed for graphene and 
the square, from scaling of the inverse participation ratio 



FIG. 5: The fractal dimension D2 for graphene and the square 
lattice as a function of off-diagonal disorder W. The D2 is 
obtained from the mean < IPR > and the typical exp < 
ln(/P_R) >) for each size L (by scaling In < IPR > and 
< ln( JP_R) > vs. InL). For W = the fractal dimension of 
graphene is D2 = 1 and for the squared lattice is D2 = 2. For 
strong off-diagonal disorder (large W) both D2 become zero 
(the zero mode state localizes). 



IPR (averaged < IPR > and typical exp < \n{IPR) >) 
vs. the linear system size L. The multifractal zero mode 
is influenced by the system boundaries even for infinite 
system size, e.g. for graphene the boundaries are respon- 
sible for the edge state structure of the Dirac point in the 
presence of off-diagonal disorder [17 1. 

A chiral system allows the presence of a mobility edge 
in 2D while the scaling theory of localization[6] forbids 
critical states in 2D orthogonal non-chiral systems. Off- 
diagonal disorder which respects chiral symmetry allows 
a topological critical state which is a zero mode in 2D. 
The problem of zero modes, apart from its intrinsic in- 
terest, is significant in other areas of physics. The first 
examples were found in the eighties and it is progressively 
realised that they are important in disordered systems, 
e.g. they can describe the low-energy physics of disor- 
dered spin chains, which can be mapped into disordered 
fermions by Jordan- Wigner transformation. Our find- 
ings support the idea that the origin of the zero modes 
is geometric. Ther multifractality is shown in two chi- 
ral disordered 2D systems, for the honeycomb and the 
square lattice. The obtained zero modes are found to be- 
have rather similarly, unless finite size effects are present. 
Their difference depends strongly on the topology of the 
perimeter, e.g in graphene for no disorder the zz edges 
on the perimeter can lead to zero energy modes. 



In conclusion, for an odd number of lattice sites a topo- 
logical zero mode exists for chiral disordered systems 
which is a mobility edge of an Anderson transition in 



2D. We have compared its multifractal properties for 
graphene and the square lattice with n.n. off-diagonal 
disorder. In both cases the isolated critical zero mode 
emerges within the localized 2D states and has ampli- 
tude only on one of the interconnected sublattices Aor B. 
The complex quantum interference in disordered systems 
which usually leads to Anderson localization is absent in 
chiral systems. In ref. 12| a perturbation theory found 
no traces of localization for chiral systems, however, non- 
perturbative effects [10| showed topological localization. 
We confirm these results. The zero mode has no topo- 
logical protection against localization which occurs for 
strong off-diagonal disorder (Z?2 —> 0). Our conclusions 
should remain valid for other chiral systems, e.g. if the 
time-reversal invariance is broken by a magnetic field and 
the invariance over spin-rotation by spin-orbit coupling. 
An old idea, whether physics can be reduced to geome- 
try, is supported by the presence of lattice topology and 
multifractality in disordered systems. In 2D topological 
phenomena (non-perturbative) exist, e.g. the quantum 
Hall efFect[l5|, etc. We have shown Anderson localiza- 
tion in chiral disordered systems appears in a topological 
sense, e.g. at the zz edges of graphene. The presence of 
zero modes is topologically protected by chiral symmetry 
but their criticality is not protected since they localise for 
strong disorder. The zero modes are also connected to 
fractional charges and monopoles[20|. The obvious ex- 
tensions of this work is to obtain the whole spectrum of 
multifractal dimensions Dq, q G [— cxd,cxd], which they 
are expected to be a non-trivial function of q, and gen- 
eralize this study to the other chiral universality classes, 

Am, ciiJUlii. 
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